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A non-Markovian master equation for a molecule interacting with a heat bath is 
obtained using the cumulant expansion method. The equation is applied to the prob- 
lem of molecular wave packet relaxtional dynamics. An exact solution is derived that 
demonstrates classical type of squeezing for the wave packet evolution and a heat-bath 
dependent frequency shift. 

1 Introduction and model Hamiltonian 

The problem of relaxational dynamics of wave packets in molecular systems is one 
of the up-to-date physical problems supported by femtosecond-time-scale spectro- 
scopic facilitiesEl as well as by single molecule sp ectroscopy.il Numerous theoretical 
investigations of the problem appeal to different approaches.il Due to the complexity 
of molecular relaxational dynamics a derivation of an exactly solving model for the 
problem has not been obtained yet. In this paper we investigate one of the models 
based on a master equation approach. 

We will start from the general potential for a single molecule interacting with 
a bath of harmonic oscillators. After the derivation of the general master equation 
we will proceed to the specific example of harmonic potential for the model, which 
admits an exact solution. 

The molecule interacting with a heat bath is separated into relevant system of 
diabatic vibronic levels and the environment of the heat bath modes £ with 
frequency u^, and creation operator b^. It is modelled by the Hamiltonian H = 
H s + H E + H SE , with 

rS 



H E = Effo;* + 1/2), 



where d^ u = \fj) (u\ is the transition operator of the system; = E^LA * s 
annihilation operator of the bath including the matrix elements of interaction 
function K^. This model was discussed previously^ and an equation of motion for 
the reduced density matrix a was derived taking the system-environment coupling 
into account by perturbation theory. In our approach to the reduced density matrix 
equation we use the cumulant expansion method.! The applicability of this method 



is rather wide because it does not appeals to perturbation theories motivation but 
relates to the statistical aspects of the influences of the heat bath. 



2 Master equation and its reduced form for the harmonic potential 

The second order cumulant expansion! gives! 

a = -^[H s ,a} + e W {-iH s t)Kexp{iH s t)a, 
Ka = H) 2 tidT([H E S(t),[H ES (T),*}}) 



where angle brackets mean averaging over environment and H ES is the Hamiltonian 
in interaction picture. This approach allows us to describe non-Markovian processes, 
when such factors as expiukx (r — t) (where uJk\ = H^ x /Ti) contain memory effects. 



We obtain for the matrix elements 

i 

J' 
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-({rxu(t)r+(r)) + ^(^(r^e^r-t)^ 
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+((rnx(r)rt (t)) + (r+ (r)r Xu (t)))a^e^ (-*)] 



Using the Markov approximation the master equation 



readsi 



Hi a 



-/C^/Cj^ {(1+n^) 5 (u^-o^) +n € <5 (u^+o^)} a KA 
+/C^ A /C Av {(l+7if) 5 (cj fcA -u; € ) (^ fcA +u; 5 )} o^] 



In the following we apply this approach to a system with a harmonic oscillator 
potential. In the interacrion picture the respective Hamiltonian of interaction reads 



H ES (t) = ^Kt(b+(t) + bs(t))(a+(t) + a(t)). (3) 

We also suggested , that there are no negative frequencies in bath. Returning to the 
Schrodinger picture we can finally write the master equation for a damped harmonic 
oscillator i 

a = —iuj [a + a, a] ([(a + + a) , a a] + [a + a, (a + + a)}) , 

+ 1) ([(a+ + a) , aa+] + [aa, (a+ + a)]) , ^ ' 

with damping rate = tc p(u)n(u>), p means density of bath states. It 

should be noted that this master equation differs from the usual form of the master 
equation for the damped harmonic oscillator!. The differences are the result of the 



full form of the molecule-bath interaction. Usually rotating wave approximation 
form of this Hamiltonian is used. The presence of the counter-rotating terms! in the 
Hamiltonian leads to the phase- dependent terms in the master equation which con- 
nect diagonal matrix elements a nn and non-diagonal ones <J mn . The later describes 
the phase-dependent effects. 



3 Method of solution and analytical results 

The solution is based on an application of the characteristic function F (A, X*,t) = 
Sp (fcr), where / = e Xa+ e~ x * a . Evaluating as usually commutators like [a, exp (Aa + )] 
one finds!: Sp([a+a, a) f) = (\*d\* — Xd\) F , Sp([a+,aa]f) = X (d\* — X) F , 
Sp([a + , aa + ] f) = —X*d\F, and etc. This allows us to transfer the operator equation 
(4) into the c-number differential equation with partial derivatives 

F = (- (tuX* + 7 (A* + A)) <9 A * + (iuX - 7 (A* + A)) d x - in {X + A*} 2 ) F. (5) 

The solution of the equation is obvious. We solve it by expanding F(X,X*,t) = 
ex P (j2m,n K mn (t) X m (— A*) n ) . Substitution F in this form into @ gives a set 
of independent systems of differential equations (SODE)i for functions K mn with 
m + n = N being a fixed number. The first two systems are: 

K 02 = -7 (n - if n) - 2 (to; + 7) K 2- 

For a wide class of initial states (coherent, thermal, squeezed, etc.) the consideration 
can be restricted by these two systems, because all higher-order elements of K mn 
are zero. Below we choose the coherent state as initial one. For this initial state the 
solutions of these SODEs read 



K w + Koi = Qoe 7 * [(j/u)) smcut + cosut] , 
if 11 + if 20 + ifo2 = n — 7 [1 + (7/w) 2 (l - cos2cjt) + (7/tD)sin2d)t] 



(6) 



where uj = \fuj 2 — 7 2 , n = (exp (hu/k B T) — 1) 1 and Qo is initial coherent state 
coordinate. The distribution of interest in coordinate space is P (Q,t)= 

I JZo dXe-^x (A, t), where X (A, t) = Sp (e iX ( a++a )a(t)) =e~^ 2 F (iX, -iX, t). In- 
tegration yields 

p(n A 2 / (Q - if 10 - if 01 ) 2 ) 

P{Q,t) = — - exp <^ - x — — — ——)■ (7) 
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I + if 20 + ifo2 + ifn) ' [ - 2 +K 20 + K 02 + K u 



Fig. 1. The plots show wave packet dynamics 7 = 0.1a;, fc^T = 3Hlj, cot = (0 : 20), Qo = 4 
(a); and the variance of the wave packet (b). 

The obtained exact solution P(Q,t) shows some kind of classical squeezing 
because the width of the wave packet slightly oscillates in time instead of it's mono- 
tone increasing for usual damping processes. This squeezing which is displayed in 
Fig.l has not been derived previously. It should be stressed that this is not quantum 
squeezing because this width is never smaller than the ground state width. Eq.6 
demonstrates a decrease of the effective harmonic oscillator frequency due to the 
phase-dependent interaction with the bath. This prediction is also new. 

The consideration of more complex potentials on the basis of proposed general- 
ized master equations is now in progress. 
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Fig. 1. The plots show wave packet dynamics j = O.lw, ksT = 3fiw, wt = (0 : 20), Qq — 4 



